Abstract We determine the finite size corrections to the large deviation function 
of the activity in a kinetically constrained model (the Fredrickson- Andersen model 
in one dimension), in the regime of dynamical phase coexistence. Numerical results 
agree with an effective model where the boundary between active and inactive 
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regions is described by a Brownian interface. 
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1 Introduction 

Glassy phenomena have proven difficult to understand: they present a variety of 
features - slow dynamics, ageing, dynamical heterogeneity, frustration - which 
make their study arduous from a theoretical point of view (see [1] for a recent 
review). Kinetically Constrained Models (KCMs) are a simple class of lattice gases 
whose dynamics shares features similar to those of glassy phenomena, with the 
advantage that no disorder is present in the model - which makes them easier to 
study (see [2,3] for reviews on KCMs). 

There is a variety of KCMs (see section 2.1 for a concrete example) which all 
share a common feature: their static properties are trivial and their complexity 
(like the phase transition phenomena) is hidden in their dynamical behaviour. This 
raises the problem of finding relevant physical parameters in order to describe and 
classify the properties of these models. As the glassy systems are characterised 
by a mixture of frozen and mobile areas, the "activity" of the system (namely the 
number of local updates during a time interval) has been proposed as a relevant 
dynamical parameter and a dynamical approach has been recently devised in order 
to define a suitable notion of dynamical free energy [4,5,6,7,8]. In this dynamical 
framework, the role of the free energy is played by the large deviation function of 
the activity. 

1.1 The large deviation function and its singularities 

For lattice gases (or more generically for Markov process with discrete configura- 
tion space) with continuous-time dynamics, the simplest definition of the activity 
is the number of configuration changes presented by an history of duration t [4, 
9] (see [10,11,12] for alternative definitions in systems with continuous degree of 
freedom) . For each history of duration t of the system, the activity will be denoted 
by the observable K t . KCMs and other glassy systems present "dynamical hetero- 
geneities" [13], i.e. regions which remain frozen during a long time. This feature 
can be captured by the probability distribution function of the activity as some 
histories present slow or inactive intervals with higher probability than in non- 
glassy systems. In the large-time limit, the probability of observing an atypical 
value K t = kt of the activity scales as 



In the infinite size limit, the function ir(k) may no longer be analytic. This can be 
interpreted as a signature of the dynamical heterogeneities [4, 14] . We postpone a 
more quantitative discussion of the singularities to section 4.3. 

From a practical point of view, it proves in fact easier to make a Laplace 
transform and consider instead [5] 



where the average is taken over histories of duration t. The parameter s plays a 
role similar to the inverse temperature in the canonical ensemble of equilibrium 
statistical mechanics: fixing s boils down to fixing the average value of the activity, 
in the same way as fixing the temperature determines the average energy. The 
functions i(>(s) and ir(k) are related by a Legendre transform: ip(s) — sup fc {7r(fc) — 
sk}. 

It has been shown for several KCMs that ip(s) presents a singularity at s J, 
in the infinite size limit, which corresponds to a dynamical phase coexistence 



Prob[if t = kt] 
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between active histories (at s < 0) and inactive histories (at s > 0) [6,7], in 
the same way as singularities of the free energy correspond to transition phases in 
equilibrium statistical mechanics. Similar singularities have been observed in other 
glassy systems (see [10,12] for binary Lennard- Jones mixtures), but the question 
of finding a generic relation between glassy properties which hold at s = and the 
singularity of t(>(s) is still open. Indeed, for any finite size system the cumulants 
of the activity: j(K™) c = (— l) n d™tp(s)\ a =o can be obtained from the function 
ip(s), but this correspondence does not hold in the infinite size limit and one may 
wonder if the singularity of ip at s = has an impact for the physics of finite 
size dynamics. Thus it is a natural question to understand how this singularity 
is build up when the system size diverges and in this article we are interested in 
the finite-size scaling of the large deviation function ip(s), especially around the 
transition. 



1.2 Finite size scaling of the large deviation function 

It has been shown that finite-size effects capture non-trivial physical features of the 
typical configurations of the system giving rise to the atypical deviation, such as 
the stationarity or the stability of the density profile in one-dimensional transport 
systems [15,16,17] or the geometrical features of the system configurations in 
glassy systems [8,12]. From a broader point of view, other quantities than the 
activity (such as the time and space integrated current) present a large deviation 
function which becomes singular in the infinite size limit [16, 18, 19], also describing 
a dynamical phase transition. 



2 Model and description of the coexistence of active and inactive 



2.1 Large deviations of the activity in the Fredrickson- Andersen model 

We focus on a one-dimensional version of the Fredrickson- Andersen [20] model (FA 
model), in periodic boundary conditions. It consists of a lattice of size L described 
by occupation numbers n = (nj)i<j<£ with sites and L identified. Each site i 
is either occupied (or 'active', rii = 1) or empty (or 'inactive', rii = 0). Transition 
rates are 



with Ci = rii-i + rii + i. The kinetic constraint Ci encodes the "dynamical fa- 
cilitation" rule: active regions favour activity in their vicinity. Compared to the 
unconstrained system (Ci — 1), the kinetic constraint does not modify the steady 
equilibrium state: each site has a Poissonian occupation number of density c, ex- 
cepted that the configuration where all sites are inactive is dismissed. The kinetic 
constraint however modifies the dynamical relaxation of correlation functions [20] 
with features similar to those of glassy systems. 

It has been shown that the dynamical free energy i/jl(s) of a system of size L, 
defined as 



regions 




(1) 
(2) 



(3) 
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presents a first order transition in the infinite size limit [6,7]: 

1 . . , f> if s < 

L i-s-oo 10 if s > 

Note that the mean activity j(K t ) s — '^t^ik^t , is also j{K t ) s = —ip' L (s). The 
transition can be interpreted as follows: 

— s < corresponds to histories where the mean activity (K t ) s is larger than the 
typical one. For these histories, the number of active sites remains extensive 
with the system size. 

— s > corresponds to histories where the mean activity {K t ) s is smaller than 
the typical one. For these histories, the number of active sites becomes finite in 
the large size limit. In particular for infinite s the value of the large deviation 
function is given by (the opposite of) the escape rate with only one active site 
in the system [7] : 

lim lim 4>l(s) — —r^ = —2c (5) 

s — ^oo L — ^oo 



2.2 Finite-size scaling and an interface model 



T 



t 



X-(j) 








Fig. 1: Model for the space-time configuration of the system in the interfacial 
regime A > A c . An island of activity density K = 4c 2 (1 — c) is delimited by two 
non-crossing biased random walks x+(t) and ar_(r), constrained to start at and 
end at at time t. 



In [8] two of the authors have considered a different scaling regime by focusing 
on values of s of the order s = X/L and for the function 

VL {X) = ^ L {X/L) (6) 

they conjectured that there exists a critical value A c > of A such that 

V (X)= lim^(A) = {-^ (7) 
L^ca —2j for A > A r 
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where K = j-^(K t ) = 4c 2 (1 — c) is the mean average activity in the system, and £ 
is the surface tension accounting for the cost of maintaining an interface between 
an active and an inactive region in the system for a long time. In [8], the limit (7) 
was derived only in a range of values A < Ao and A > Ai for some parameters 
< Ao < Ai and not up to the conjectured critical value A c = J?. 

Note that S ^ r^: the typical configurations of the system at finite A > A c 
are not given by those of the s — > oo limit. In particular, they present more than 
a finite number of active sites. Our aim in this article is to identify the typical 
configurations occurring at A > A c , and to determine the finite size corrections 
that they imply on the infinite size result (7). These configurations are interesting 
to characterise because they are the first to appear when increasing A (that is, 
they are the first to appear when considering histories of the system displaying an 
activity K t lower than the typical one). 

We introduce now a simplified dynamics in order to model the configurations at 
A > A c . In the slow activity regime, the system can be described at a macroscopic 
level by a small active "island" of mean activity K in a large sea of an inactive region 
(see figure 1). In the unbiased dynamics (A = 0), the inactive region would be 
invaded and become active. Thus at the macroscopic level, an interface between an 
active and an inactive region should perform a biased random walk with effective 
jump rates p, q which take into account the growth of the active region. When 
A > A c , the growth of the active region is penalised as the activity of the system 
is proportional to the area of the active droplet. 

More precisely, the boundaries x + (t) and of the active region perform 

non-crossing random walks of jump rate p (resp. q) to the left (resp. right) for x + (t) 
and mirror rates for x~(t). For simplicity the walks are constrained to start from 
x = at time 0, and to come back to at final time t (this assumption does not 
change the large time asymptotics). In this effective description, the total activity 
in the system is proportional to the area of the active droplet and approximated 
by K / Q dr [x + (r) - x_ (r)] with K = 4c 2 (1 - c) the mean density of activity. Thus 
the counterpart of (e~ sKt } reads 



This is the simplest model one can think of to represent the separation between 
active and inactive regions in the system. In particular adding more interfaces 
would lead to a metastable situation where the active regions eventually merge 
together to form a unique island of activity. We think that the interface model 
represents the correct dynamics of the system at large scale, but we have not 
found a rigorous derivation starting from the microscopic dynamics. However, 
the numerical results of section 3 support the scaling derived from the simplified 
model (11). 

Thus we conjecture that the finite size corrections to the large deviation func- 
tion (/ 3 l(A) (6) for A > A c are related to 






where (•) denotes the average over trajectories x±(r)o< r <t without constraint 
at final time. 



(fi L (X) = lim i log Z eS {±t) 



(9) 



Inspired by the study of interfaces in the static Ising model [21,22], and using 
results from Brownian bridge theory [23], we show in appendix A that this leads 
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Fig. 2: Evaluation of the large deviation function </?l(A) using the cloning algo- 
rithm (blue circles, increasing sizes L £ {8, 16, 32, 64} from bottom to top at pos- 
itive A), and using direct diagonalisation of the operator of evolution (54) (plain 
green line, L — 8 run as a check). The red dashed line is the infinite L result — KA 
for A < A c , while the purple dotted horizontal line is the infinite L result — E for 
A > A c . We took c = \. 



to the following scaling at large L 



AIK 

<p L {\) = -4^q( — — ) 2-tai (10) 



where a% ~ 2.3381... is the first zero of the Airy function on the negative real axis. 
As a consequence, we expect that the finite size scaling of the microscopic model 
should be given by (10) plus the extra cost — £ for creating the interfaces 

2 

^ L {X) = -S-A^r q (-^-\ 2-iai (11) 



for appropriate choice of the effective parameters p, q (see section 4.1 for a dis- 
cussion on the effective jump rates). In other words the interface model we have 
considered leads to L~i corrections to the constant —S. 



3 Numerical results 

3.1 Results from the cloning algorithm (1): the free energy 

To investigate whether the finite size corrections (11) inferred from the inter- 
face model are correct, we have measured v'l(A) in numerical simulations. Since 
large deviations are by definition difficult to measure, a direct sampling of iPl(s) 
through (3) is not achievable. We have resorted to a continuous-time version [24, 
25] of the Giardina-Kurchan-Peliti cloning algorithm [26] in which the dynam- 
ics is modified so as to make the large deviation typical, at the price of muta- 
tion/selection rules between a large number of copies of the system (see [27] for 
a review on cloning algorithms). Those algorithms have already been used to de- 
termine large deviation functions in lattice gases [6,28,29] but not in the scaling 
regime s — X/L that we consider in this article. 

A first result (figure 2) is that the large deviation function agrees qualitatively 
with the conjectured infinite size result (7): the large deviation function ^l(A) 
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Fig. 3: (Left) Log-log plot of (Pl(\q)+Z!, for S — 0.077, at fixed Ao = 4.6 as a func- 
tion of L: the numerical evaluation (blue dots) fits with a power law corresponding 
to the exponent exponent a = | (red line). (Right) Plot of L Q (</?i(A) + S) for 
different values of L (L e {64, 96, 128, 160, 192, 256, 320}). The curves collapse on 
a single master curve — 1.05A" (dashed blue), in agreement with the interfacial 
model result (11). The parameter c of the model is c = |. 



tends to become linear for A < A c and constant for A > A c as L increases. The 
critical value is determined as A c = UfK. 

The scaling of the deviations from the infinite size result is examined in figure 3. 
In agreement with the interfacial model result (11), ¥>l(A) + £ scales in L~ a at 
fixed A (figure 3, left), with a = |, while the rescaled curves L a (ipL(X) + E) 
collapse onto a master curve — AiA Q , with A\ ~ 1.04 (figure 3, right). Simulations 
were performed at mean density c = |, but the results and the scaling analysis we 
present do not depend on this value. 



3.2 Results from the cloning algorithm (2): the density 

One may test another consequence of the interfacial model by computing the 
density of active sites. Taking the derivative with respect to s in (3) leads to 
\{K t ) 8 = — tP'l( s )- Thus using the relation (6) implies that ~{K t )\ = —L(p' L (X) 
and from (11) the mean activity j{K t )\ for histories weighted by e~ XKt ^ L scales 
as Li. Therefore, one expects that the average width of the active droplet is 
Ls when A > A c . Strong finite-size effects are still present (figure 4) and wouldn't 
allow to check precisely the power | of the scaling relation (11). To understand the 
origin of these corrections a useful tool is the escape rate r(C) from a configuration 
C = («i)i<i<L defined as the sum of the jump rates from C: r(C) = Ylf=i [(1 ~ 
c)rii + c(l — m)] (rii-i + rii + i). As shown in appendix B, the fluctuations of K t and 

of the time integral of the escape rate R t = J Q dr r(C(r)) are closely related (this 
result is valid in general): 

i> L (s)= 1 t (Kt)s-~{R t )s (12) 

In the interfacial regime A > A c , <^l(A) goes to a constant — E when L goes 
to infinity, while both j(K t )\ and j(Rt)x grow with L so that (12) describes the 
cancellation between those growths. More precisely, assuming the scaling form 

</?l(A) = —E + L^ a Lp 1 (\) + o(L~ a ) (13) 
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Fig. 4: In red (resp. green): log-log plot of j(K t )\ (resp. j(Rt)x ) a t fixed A = 4.6 
as a function of L (same values as in figure 3). Dots are numerical data while lines 
are power law asymptotics with exponent 1 — a = |. Duplicated dots correspond 
to different runs. 



one sees by differentiating (12) that both (K t )\ and (Rt)x scale in the same way 
with L 

~<2r t ) A ~L 1 -«* 1 (A) + o(L 1 - Q ) i( J R t ) A ^L 1 - Q r 1 (A)+ (i 1 -«) (14) 

where the exponent 1 — a comes from the relation ^(K t )\ = —Lip' L (X) and r\ =k\. 
We also note that at c = i, r(C) is (twice) the total number of active sites 
so that the \{Rt)x also represents the mean density of active sites at fixed A, 
which thus scales as L x ~ a = L$ (we expect the result to hold also for c / ^ 
although -j(Rt)x is not the mean density anymore). We conclude that a numerical 
check or evaluation of the exponent a is better done on </5_l(A) than on the mean 
activity, the density or the escape rate, since those quantities present large finite- 
size corrections e.g. of the form 

~(K t ) x ~ L x - a MA) + L'<k 2 {X) + C K + 0(L- a ) (15) 

~(Rt) x ~ L 1 ~ a n(A) + L^r 2 (A) + C R + 0(L- a ) (16) 

where by cancellation from (12) ki — r%, k 2 — r 2 and Ck — Cr = 



3.3 Results in an independent site approximation 

We now consider another approach to compute the dynamical free energy (^(A). 
It can be shown that ^l(A) is the largest eigenvalue of a symmetric operator of 
evolution W A ym [5] (see also appendices B and C) which acts on the vector space of 
all configurations {|n)} of the system. We thus have the Courant-Fisher equality 

(x\w s r\x) 

^ = ffijg, (x\x) (17) 



where \X) — ^2 n X(n)\n) is a vector. 
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L a (<p L (\,L) + E) 



Fig. 5: Plot of the scaling function L a ((pL(A, L) + E) for different values of L 
(L € {32, 36, 40, . . . 96, 100}) determined from the Bernoulli optimisation (18) (blue 
dots). The curves collapse on a single master curve — AiX a (dashed red). The 
values of E and Ai are quite different from those of the original FA model: here, 
E w 0.251 and A x w 0.897. 



Restricting the optimisation on laws X(n) representing products of Bernoulli 
distributions for evaluating the largest eigenvalue of W| ym , one obtains in ap- 
pendix C the following estimation for the large deviation function: <Pl(X) < 
(pf ein (\) with 



BarnMX El<i<i e V^ 1 ~ C )Pi^ ~ Pi) ~ C ( X ~ Pi) ~ C 1 ~ C )Pi] + ft-l) 

Vl (A) = max — = r 

In} 1 ~ IL^KiA 1 ~ W 

(18) 

where < ^fpl < 1 is the parameter of the Bernoulli law at site i. 

The mean value of the mean fraction of active sites L dr Jii(r) for his- 
tories weighted by e~ sKt is given by 



L) = ^(e-^i ^r^n,(r))/(e-^) = (i|i^n,|i?) (19) 



where \L) and are the left and right eigenvectors of W s associated to the 
eigenvalue </?l(s). The corresponding value in the Bernoulli projection reads 



p Bern (A,£) 



1 



1 



L l-Ui<i<L^-pi) , 



<i<L 



(20) 



for a set {pi} which optimises (18). Numerically solving (18) one obtains p Bern (A, L) 
and the large deviation function. The collapse for (p Bern (A) works again with a = | 
(see figure 5), while again p Bern (A, L) displays stronger finite-size effects. 

We expect the true minimiser to be very different from the approximation 
by a Bernoulli distribution. Indeed the sites located at the interface should be 
very correlated. This explains why the values of E and A\ measured within the 
Bernoulli approximation are different from those of the original FA model (see 
section 4.2). 
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4 Discussion 

4.1 Determination of the effective jump rates p and q, and the surface tension E 

The exponent a = | observed numerically matches the one predicted by the 
interface model (11). One can also compare the value of the coefficient A\ in the 
scaling 

^ L (\) = -£-A 1 (^y + o(L~ a ) (21) 

In a very crude approximation, where correlations are neglected, one can imagine 
the interface as a single active site at position x(t) separating a region with only 
inactive sites from an active region sampled according to a Bernoulli measure 
with parameter c. It this case, the parameters p, q of the interface model defined 
in section 2.2 can be estimated as follows. We have q = c: the interface grows at 
rate c by activating a site on the left of x(t). On the other hand, the value of p 
may be estimated to p — c(l — c): due to the kinetic constraint the first active site 
is inactivated with rate 1 — c provided that the site to its right is active, which 
occurs with probability c. From (11) and K = 4c 2 (1 — c), this yields 

2 

A 1 =4Jp^(^^-) 2-3ai=K3ai (22) 

and thus A x ~ 1.312 for c = |. This is close to the value A x w 1.05 ± 0.01 
observed numerically, the discrepancy arising in part from correlations between 
neighbouring sites around the interface. 

The value of the dynamical surface tension E can be compared to analytical 
predictions. In [8] an expression of E was derived for the East and the FA models 
with fixed boundary condition. A similar expression holds for the periodic FA 
model we are interested in: 

E = - lim sup(\/p|W sym (i)|\/p) (23) 

where the supremum is performed over all probability distributions on the set of 
configurations {n} of the system, with the condition that site 1 is inactive: m = 0. 
The vector | \fP\ is the vector of components yj -P(n). The symmetrized operator 
of evolution W sym (L) for a system of size L is defined in appendix C. Note that 
the formulation of this extremalisation principle in [8] involves a Dirichlet form 
which is equal to the expression maximised in (23), up to a sign. Defining now the 
projector Px, onto the configurations with site 1 inactive, we can replace W sym (L) 
by P L W sym (L)P L in (23) and relax the condition on P. This shows that: 

E = Urn E L with Z" L = - maxSp (P L W sym (r)P L ) (24) 

where Sp denotes the spectrum of an operator. We thus have expressed the dy- 
namical surface tension using the maximum eigenvalue of an operator, in a similar 
way as for the free energy iPl{s) in (17). The value of El was computed for sys- 
tems sizes L < 15 by direct diagonalisation (figure 6). By fitting the results using 
a reasonable form of the finite-size corrections we obtain the value E = 0.0771, 
which is numerically compatible to the one obtained from the cloning algorithm 
E = 0.077 ± 0.0005. Note that the finite-size surface tension El displays strong 
finite-size effects fitted through the form El = E + AL~z , but we have no theo- 
retical justification for the power |. 
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Fig. 6: Plot of the surface tension El at size L obtained by diagonalisation of 
the operator appearing in (24) (blue dots). To determine the surface tension E — 
liniL^oo El, we performed the fit El = E + AL~i (green line). The exponent | 
was chosen as the most likely one in regards of the numerical results. The result 
E = 0.0771 (horizontal red dashed line) is very close from the one obtained with 
the cloning algorithm (0.077, see figure 3). 

4.2 Universality 

The Bernoulli approximation developed in part 3.3 presents the same scaling ex- 
ponent a = | with different constants E and Ax (see figure 5). This tells that, 
on the one hand, the independent site approximation of the Bernoulli optimisa- 
tion (18) is quite far from being correct: the large deviations do not match at 
all. This fact was also checked by diagonalising the operator of evolution W s and 
finding the eigenvector associated to the maximal eigenvalue tpL(s), for small sys- 
tem sizes (L < 15). We found that the corresponding state can't be factorised on 
independent sites and presents correlations. On the other hand, the robustness of 
the exponent | is an indication that other models in the same class also present 
(an) interface(s) in the coexistence regime. 

Besides, let us note that the same expression as (18) appears when considering 
the field theory associated to the operator of evolution W^ ym [5] and computing 
¥>l(A) assuming a time- independent saddle-point. This tells that such an assump- 
tion is not valid. 



4.3 Finite size scaling of the large deviation function 7Tl(£;) 

At fixed system size L, the large deviation function ttl defined as 

Prob^ = Lkt] ~ e t7TL(k) (25) 

is the Legendre transform of i[>l (3) 

TT L (k) = inf {sLk + tPl(s)} , iPl{s) = swp{-sLk + 7r L (/c)} . 

s ft 

Using the parameter A = sL, the large deviation function becomes, in the large L 
limit, for reduced activities 

Vfe € [0,K], n(k) = lim n L (k) = inf {Afc + <p(\)} = -E(l - (26) 

L-s-oo A>0 ' K 
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Fig. 7: Schematic plot of the finite-size scaling of the large deviation function 
T^hik) = lim^oo logProb [K t — tLk] , deduced from the finite-size scaling of ^^(X). 
The regime of smallest k (k <C L^ 1 ) corresponds to configurations with finite 
number of active sites. The intermediate regime (L~ 3 < k <C K— L^ 1 ) corresponds 
to the regime of phase transition, where the large deviation function is linear 
in the infinite size limit (cf. equation (40)). The picture for an unconstrained 
dynamics (C\ — 1 in (1-2)) has a completely different large deviation scaling 
7r™ c (fc) = L7r" nc (A:) with inactive histories exponentially much less likely. 



where ip is defined in (7). This follows from a simple computation (see [30] for a 
review on ldf in physics and mathematics). Note that negative activities cannot 
be produced. The physical picture behind the linear behaviour (26) is a two step 
mechanism. An activity deviation of order kt over the time interval [0, t] is pro- 
duced by first blocking the system in an inactive state during a time (1 — j|)t and 
then letting the system in the stationary state (with mean activity K) during a 
time 4i. This leads to an exponential cost U(l — k)t. We stress the fact that the 
switch from the inactive to the active state occurs on time scales much smaller 
than t and therefore has no impact on the large deviation function obtained in 
the large t limit. This behaviour is characteristic of a first order phase transition 
between the active and inactive regime. We examine below the consequences of 
the finite size corrections to the large deviations within both regimes. Results are 
summarised on figure 7. 

The active regime: 

For activities larger than K, the scaling is different as the constraints do not 
play a major role and the behaviour is similar to non-constrained systems, one 
has 

Vfc > K, Tr(k) = lim ^n L (k) = inf (sfc + ip(s)} , (27) 

L— ¥00 L s<0 I J 

where tp(s) — lim£_ i . 00 -^-^(s) which leads to 

Vfc > K, Prob^ ~ tL k] - e ti#(fc) (28) 
In particular n is expected to be a smooth function for k > 

Vfc>K, n(k) = --^-(k-K) 2 + 0({k-Kf) (29) 
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where K2 = j-^(Kf) c is the variance of the activity. 

The finite size scaling asserts that the transition for iPl(s) takes place away 
from s — and therefore, we expect that for very small deviations of the activity 
below K, Ltt approximates the large deviation function 7r^(fc) = inf^o {Xk + ip^iX)}. 
This would hold only for s > i.e. k — K w ~ (by optimising in (27)). This 
means that in a small window around K of order -t, the large deviations are of 
the form (29) 

^dk) = -^{k-K) 2 + 0{L^) (30) 

The inactive regime: 

In the inactive regime, the finite size scaling takes into account the deviations 
of the width of the active droplet which will contribute to the activity deviations 
when k is close to 0. The finite size Legendre transform (26) reads 

Vfce[0,K], 7T£,(fc) = inf {Afe + ^i(A)} (31) 

where the optimal A is given by k = —tp' L (A). The finite size scaling (11) is valid 
for A > A c 

V L W = -Z-A 1 ^ 2/3 =* fc =3Z^ A ~ 1/3 (32) 
This allows to estimate 7Ti(fc) for k w and A 1 

L3 

4/13 1 

*l(*)*-B-^JW (33) 

This scaling should remain correct if the droplet width is very large microscop- 
ically, i.e. for k 1/L. For k « L~z\~3 with A close to zero, the probability 
of observing an active droplet of width kL should vanish exponentially fast with 
a rate given by (33). But we will see below that these larger droplets do not 
contribute to n^(k). 

The intermediate regime: 

We consider now the intermediate deviations in k. The two-step mechanism 
described earlier has to be slightly modified to take into account the finite size 
corrections found in the active and inactive regimes (30), (33). Thus one expects 
that 

Prob[X t ~ tL k] ^ ~ (34) 

1 



sup < ( 

>=1 i">3 



sup <j (\xp j -tXi - - tx 2 



27 L 2 kl 



where the supremum is taken over the activities k\, k 2 and xi, x 2 such that 

l=Xi+X2, k = xiki + x 2 k 2 (35) 

Suppose that k <C K — i. As k\ — K is at most of order 1/L, one can neglect the 
deviations with respect to k\ and set k\ — K. This leads to optimise over fc 2 with 

IK k 

fci=K, xi = l-x 2) a; 2 = — — €[0,1] (36) 
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Thus (34) reads 

Probst ~ tL k] ~ exp ( -i(K - k) inf {G(fc 2 )} ] (37) 

t— >oo \ k 2 <k J 

where 

1 4A 3 1 

The function G reaches its minimum for 

^A 2 L e E 3 (A + K 2 L 2 i7) + AKL 4 E 2 ^j — AL 2 i7 /2A\ 1/3 1 



L 2 EVy/A 2 L e E 3 (A + K 2 L 2 E) + AKL 4 E 2 



4:A' 

with A — Finally we get 



3 



E J L 2 / 3 
(39) 



r o urrz *r n j — (K — fc) G(k*), for ^<fc«K-i 

hm - log Prob ifj ~ tL k ] = < , : , ' , L (40) 

i^oot J 1 -(K-fe) G(fc), for y ' 

Thus for a deviation k > k*, the active droplet (in the inactive phase) has a 
microscopic width located around Lk* and the large deviation is still linear with 
minor corrections compared to the limiting case. For k < k* then the system 
remains all the time in the inactive phase (x% = 1) and the droplet width shrinks 
leading to a non-linear large deviation cost in k. The fact that larger droplet 
widths cannot be observed in the large deviations is due to the first order phase 
transition. 

The very inactive regime: 

For s — > oo (i.e. A ^> L) a single remaining site is active and one obtains the 
following asymptotics [6]: 

iPl(s) = -Too + 2ce~ s + 0{e~ s ) for s -> oo (i.e. A > L) (41) 

where = 2c is the mean escape rate in the configuration with one active site. 
Performing the inverse Legendre transform, this implies that 

n L (k) = - roo +Lk(l-log^) forfc«L- 1 (42) 

2c 



5 Conclusion 

We have shown that the dynamical phase coexistence occurring in the Fredrickson- 
Andersen model for histories at small positive s is well described by two Brownian 
interfaces enclosing an island (or "space-time bubble" [13]) of activity of width 
in a system of size L. The scaling of this physical picture is reflected in the finite 
size scaling of the dynamical free energy iPl(s) of the model. We expect the same 
picture to be valid for a wider class of one-dimensional kinetically constrained 
models where the particle number is not conserved. It would be interesting to in- 
vestigate the relation between our result and the statistics of large inactive bubbles 
in the non-modified (s = 0) dynamics. 

In general, the finite size scaling exponents depend on the dimension and on 
the nature of the constraints [8] . Thus it would be interesting to extend our study 
to more general dynamics, in particular, to understand the interface fluctuations 
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in higher dimensional models. This would be key to connect our approach to 
realistic glass formers. The quantitative link between those realistic (in general, 
atomistic) models and kinetically constrained models have been explored in a va- 
riety of studies [31,32,33,34]. Although it is beyond the scope of this article to use 
this correspondence to provide quantitative predictions on activity large deviation 
function in realistic models, we expect that the finite-size scaling exponents are 
related to geometrical features of inactive regions (see [12] for an example of such 
exponent). 
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A Scaling of i/)l(s) m the interface model 

We determine in this appendix the finite size corrections to the function <^_l(A), defined 
in (9), in the effective interface description discussed in part 2.2. Let us first focus on 
a system with one boundary x(t) between the active and inactive regions. It performs a 
random walk of jump rate p (resp. q) to the left (resp. right), starting from x = at time 
0, and constrained to come back to at final time t. The walk x{t) is reflected at so that 
x(t) > 0. Appropriate values of p and q are discussed in section 4.1. The computation is 
done for generic values of p and q. 

In the original model, the cost of maintaining the interface (that is, for x(t) to come 
back in at time t = t) is given by the surface tension E. In this effective description 



e 



;/ 'dr*(r) S(x(t)=0)) 



Z -» IM)= — <.(,(.)-»„ ,43) 

where K = 4c 2 (1 — c) is the mean density of activity and (•) denotes the average on 
trajectories x{r)o< T <t without constraint at final time. The normalisation (S(x(t) = 0)) p 
is e.g. fixed from Z e g(0,t) = 1. 

Let us denote by P(x, X, t) the probability of being in x at time t, having observed 
a value X of the area f dr x(t). The initial condition is P(x, X, 0) = 8 x ,oS{X). Defining 
P(x,s,t) = J dX e~ s x P(x,X, t), one has (using the constraint at final time x(t) = 0): 

= ££^ (44) 
P(0,0,t) 

Moreover, from Feynman-Kac formula, the time evolution of P(x, s, t) is given by 

d t P(x,s,t) = pP(x + l,8,t) + qP(x - l,s,t) - {p + q)P(x,s,t) - sKxP(x,s,t) (45) 

and a reflection at x — 0. To symmetrize the walk, we now set Q{x, s, t) = (q/p)z P(x, s, t). 
One has again Z e g(s,t) — SL^'f and the evolution of Q(x, s,t) writes 

dtQ(x, s,t) = ^/pq^Q(x+l, s,t) + Q(x — 1, s, t) — 2Q(x, s, t)] 

- sKxQ(x, s,t)+ (2 v / w - p - q)Q(x, s, t) (46) 



The normalisation has changed but we see that apart from the constant loss rate (2^/pq — 
P ~ q) i Q describes a symmetric walk with the same term sKxQ(x, s, t) corresponding to 
weighting trajectories by the area J dr x(t). In other words, defining at last Q(x,s,t) = 
Q{x,s,t)e- t{2 ^*- p - q) , one has 



(47) 
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and from the equation of evolution 

dtQ(x,s,t) = v / pq[Q(x + l,s,t) + Q(x- l,s,t) -2Q(x,s,t)j - sKxQ(x,s,t) (48) 

we see that Q(0,0,t) does not increase exponentially in time since for s = the equation 
describes a simple symmetric random walk. Moreover, in the large size limit (s — ¥ 0), the 
evolution of Q{x, s, t) is governed by the continuous in space operator 

- sKx = [iflg - (49) 

with reflecting boundary condition in 0. For bridges, the spectrum is known and its largest 
eigenvalue is given by [23] 

i7' PE » = 2^ 3 2-'ai (50) 

where a.\ ~ 2.3381... is the first zero of the Airy function on the negative real axis. In 
periodic boundary conditions, one has two interfaces and the equivalent jump rates are 
multiplied by 2. Finally, this gives in the s — > limit: 

2 

^(^^(^V^r (51) 



B A generic identity between large deviation functions 

In this appendix we prove an identity used in part 3.2 between large deviation functions 
associated to the activity Kt and to the escape rate Rt defined below. We consider a 
Markov process on a finite number of configurations {C}, with transition rates W(C — > C ) 
between configurations. The activity Kt is an history-dependent observable increasing by 
1 upon jumping from C to C' . The probability P(C,K,t) of being in C at time t having 
observed a value K of the observable Kt thus evolves in time through 

d t P(C, K,t) = J2 W ( C ' -> c ) p ( c '. K-l,t)- r(C)P{C, K, t) (52) 
c 

with r(C) = J^c' W(C? — > C') the escape rate from configuration C. The Laplace transform 
P{C, s, t) = Ex e- sK P{C, K, t) verifies 

d t P(C, s,t)=J2 e- a W(C' -> C)P(C', s, t) - r(C)P(C, s, t) (53) 
c 

The cumulant generating function ip(s) defined in the infinite time limit as {e~ sKt ) ~ ^i a > 
is the largest eigenvalue of the operator W s of elements 

(W 3 ) cc , = e~ a W{C -+C)- r{C)5 cc > (54) 

since (e~ sKt ) — P(C, s,t). Both Kt and R t — f drr(C(r)) quantify the activity of the 
histories. Their large deviation functions (ldf) are closely related. Indeed, let's consider the 
joint ldf 

¥(s,o)= lim \\og{e- sKt -' 7Rt ) (55) 

As previously, one checks that ^(s, <r) is given by the maximum eigenvalue of the operator 
Ws.o- of elements [7] 

(W S , CT ) CC , = e~ s W(C -)• C) - (1 + <T)r(C)S CC ' (56) 
It verifies the symmetry W S , CT = (1 + o")Wj,+i og (i+ (T ) i o an d so does the ldf: 

f (s, a) = (1 + a)<P(s + log(l + a), 0) (57) 
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Besides, the mean values of Kt and Rt in the s-state are given by 

l(K t ) s = -d s V(s,a)l =0 1 -{R t ) 3 = -d^{s,a)\^ (58) 
Differentiating the symmetry (57) with respect to a and sending a to 0, one gets 

ft^MUo = *(a,0) + d s *(s,a)l =0 (59) 

which implies 

^(s) = j(K t ) s -](Rt)s (60) 

This relation is generic. It leads to the relation (14) between scaling exponents for the FA 
model in the inactive regime s — X/L. 



C Bernoulli approximation to determine cpl(\) 

In this appendix, we obtain the expression of the Courant-Fischer optimisation principle of 
part 3.3 for Bernoulli states. Before this, one needs to symmetrize the evolution operator 
W 3 introduced in (54). We take the notation of appendix B. Assuming that the jump 
rates verify the detailed balance symmetry W(C — > C')P ecl (C) = W(C — > C)P eq (C), it 
is generically possible to symmetrize the operator of evolution W s through the similarity 

transformation W| ym = P e q 2 W s P e q, where P eq is the diagonal operator of elements P e q(C). 
Upon symmetrisation, we have that tp{ s ) is a ^ so the largest eigenvalue of the symmetric 
operator W| ym of elements 

(W s / m ) cc , = e- s [W{C -> C)W(C -> C')] * - r(C)8 cc , (61) 

It is convenient to represent the operator of evolution in terms of spin = operators 

a and n. On each site i, is the creation/annihilation operator and n is the counting 
operator. They are defined by 

cr+|0) = |l) cr~|0)=0 rl|0> = (62) 

cr+|l)=0 cr _ |l) = |0> n\l) = |1) (63) 

Where |0) and |1) are the vectors for empty and occupied states. One has 

W s = J2 [ e ~" fat + 0- ~ fai) -c(l-fk)-(l-c)hi](n i+ i+ni-i) (64) 

l<i<L 

Transition rates obey detailed balance with respect to the product Bernoulli measure of uni- 
form density c (conditioned to exclude the fully inactive configuration). The symmetrized 
operator of evolution writes 

W s / m = Yl [e" Vc(l - c)(a+ + err) - c (l - ft.) - (1 - c )ft 4 ] (ft i+1 + ft._ x ) (65) 

l<i<L 

Defining the vector \p) corresponding to the Bernoulli distribution of density ^fp: 

|p) = ^|l) + (l-^)|l) (66) 

and taking \X) in (17) to be the state |pi . . . Pl) c ° u conditioned to exclude the fully empty 
configuration: 

\pi . . . p L ) c ° nd — ^2 (ni . . . n L \pi ■ ■ ■ Pl) \ni . . . n L ) (67) 
one obtains (18) by direct computation. 
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